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An analytical method of studying strong long-range electron-phonon and Coulomb interactions in 
complex lattices is presented. The method is applied to a perovskite layer with anisotropic coupling 
of holes to the vibrations of apical atoms. Depending on the relative strength of the polaronic shift, 
Ep, and the inter-site Coulomb repulsion, Vc, the system is either a polaronic Fermi liquid, Vc > 
1.23 Ep, a bipolaronic superconductor, 1.16 -Ep < Vc < 1.23 Ep, or a charge segregated insulator, Vc < 
1.16 -Ep. In the superconducting window, the carriers are mobile bipolarons with a remarkably low 
effective mass. The model describes the key features of the underdoped superconducting cuprates. 
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There is clear experimental and theoretical |^-|l^] 
evidence for strong electron-phonon (el-ph) interaction 
in high-Tc cuprates. Electron correlations are also im- 
portant in shaping the Mott-Hubbard insulating state of 
parent undoped compounds . The theory of high-Tc 
cuprates must treat both interactions on equal footing as 
was suggested some time ago |6). In recent years many 
publications addressed the fundamental problem of com- 
peting el-ph and Coulomb interactions in the framework 
of the Holstein-Hubbard model [pl]-p5| where both inter- 
actions are short-range (on-site) . The mass of bipolaronic 
carriers in this model is very large and the critical tem- 
perature is suppressed down to a kelvin scale. However, 
in the cuprates the screening is poor so that the el-ph 
interaction necessarily has to be long range. Motivated 
by this fact, we have proposed that a long range Frohlich, 
rather than short range Holstein, interaction should be 
the adequate model for the cuprates |jr^,|l8|. A small 
polaron with the Frohlich interaction was discussed long 
time ago Analytical and exact Monte-Carlo 

[0 studies of the simple chain and plane models with a 
long-range el-ph coupling revealed a several order lower 
effective mass of this polaron than that of the small Hol- 
stein polaron. Later the polaron and bipolaron cases of 
the chain model were analyzed in more detail in Refs. 
pO| and confirming low masses of both types of car- 
riers. Qualitatively, a long-range el-ph interaction results 
in a lighter mass because the extended lattice deforma- 
tion changes gradually as the carrier tunnels through the 
lattice. 

In this Letter we study a realistic multi-polaron model 
of the copper-oxygen perovskite layer, a major structural 
unit of the HTSC compounds. The model includes the 
infinite on-site repulsion (Hubbard U term), long-range 
inter- hole Coulomb repulsion Vc, and long-range Frohlich 
interaction between in-plane holes and apical oxygens. 



Wc find that within a certain window of Vc the holes 
form intersite bipolarons with a remarkably low mass. 
The bipolarons repel and the whole system is a super- 
conductor with a high critical temperature. At large Vc, 
the system is a polaronic Fermi-liquid, and at small Vc it 
is a charge segregated insulator. 

To deal with the considerable complexity of the model 
we first describe a theoretical approach that makes the 
analysis of complex lattices simple in the strong cou- 
pling limit. The model Hamiltonian explicitly includes 
long-range electron-phonon and Coulomb interactions as 
well as kinetic and deformation energies. An implic- 
itly present infinite Hubbard term prohibits double oc- 
cupancy and removes the need to distinguish fermionic 
spins. Introducing spinless fermion, Cn, and phonon, 
dma, operators the Hamiltonian is written as 

H = -Yl [^(" - n')4cn' - Vein - n')cj,c„c|,,c„. 
- ffg (m - n) (e 

nm 

+ MLcrfma + ^j . (1) 

ma ^ ^ 

Here is the polarization vector of ath vibration co- 
ordinate at site m, Um-n = (m — n)/|m — n| is the unit 
vector in the direction from electron n to ion m, and 
ga (m — n) is a dimensionless el-ph coupling function. 
[ga{ui— n) is proportional to the force acting between 
m and n.] We assume that all the phonon modes are 
dispersionless with frequency lu and that the electrons 
do not interact with displacements of their own atoms, 
ga{0) = 0. We also use Ti = 1 throughout the paper. 

In the limit of strong el-ph interaction it is conve- 
nient to perform the Lang-Firsov canonical transforma- 
tion H . Introducing S = J2mna9a{m ~ n){e^a ■ 
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Um~n)cJiCn(<ijiiQ — dma) One obtains a transformed 
Hamiltonian without an explicit el-ph term 



w(n - n')c]jCnC^/Cn' - -Bp ^ c|jCn. 



(2) 



The last term describes the energy which polarons gain 
due to el-ph interaction. Ep is the familiar polaronic 
(Franc-Condon) level shift 



(3) 



which we assume to be independent of n. Ep is a natural 
measure of the strength of the el-ph interaction. The 
third term in Eq. (H) is the polaron-polaron interaction: 



v{n - n') = 14(n - n') - V^d^^ - n'), 



(4) 



Vpa(n - n') = 2cj^5a(m - n)ga{m - n') x 

ma 

(^ma ■ ^m— n)(^mQ ' ^m— n')' (^) 

where Vpa is the inter-polaron attraction due to joint in- 
teraction with the same vibrating atoms. Finally, the 
first term in Eq.(0) contains a transformed hopping op- 



<5-nn' = T{n - n') exp 



^ [5a (m - n)(ema ■ Um-n) 



- ga(m - n')(ema ■ Um-n')] (c'La ^ dma)] ■ (6) 

At large £'p/T(n— n') this term is a perturbation. In 
the first order of the strong coupling perturbation theory 
|§ , CTnn' should be averaged over phonons because there 
is no coupling between polarons and phonons in the un- 
perturbed Hamiltonian [the last three terms in Eq.(||)]. 
For temperatures lower than the result is 



t{n - n') = (o-^ 



Iph 



r(n-n')exp[-G^(n-n')], 



(7) 



This is the simplest way to calculate and (bi)polaron 
masses once the 'static' parameters Ep and Ipa are 
known. 

It is easy to see from the above equations that the long- 
range el-ph interaction increases Ep and I^pa but reduces 
(when measured in natural units of EpjiS). Thus 
polarons get tighter and at the same time lighter. Bipo- 
larons form when V^pa exceeds 14 and they are relatively 
light too. We note that the Holstein model is the limit- 
ing case with the highest possible G^ — Epjuj. In this 
respect, the Holstein model is not a typical el-ph model. 

To obtain an analytical description of the multi- 
polaron system we restrict our consideration to the strong 
coupling case t < \v\. In this regime the polaron kinetic 
energy is the smallest energy and thus can be treated as 
a perturbation. The system is adequately described by a 
purely polaronic model: 



G^(n- n') = ^^^(m- n)(ema • Um-n) x 

ma 

[5a(m - n)(ei„Q • Um-n) - 5a(m - n')(emQ • Um_„/)] . (8) 

By comparing Eqs.(||), (||), and (||), the mass renormal- 
ization exponents can be expressed via Ep and Vpa as 
follows 



G2(„_n') = ;^(i?p-iT/pa(n-n') 



(9) 
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-ffpcrt = - ^ t(n - n')4c„/. 
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The many-particle ground state of Ho depends on the 
sign of the polaron-polaron interaction, the carrier den- 
sity, and the lattice geometry. Here we consider a two 
dimensional lattice of ideal octahedra that can be re- 
garded as a simplified model of the copper-oxygen per- 
ovskite layer, see Figure 0. The lattice period is a = 1 
and the distance between the apical sites and the cen- 
tral plane is h = a/2 = 0.5. All in-plane atoms, both 
copper and oxygen, are static but apical oxygens are in- 
dependent three-dimensional isotropic harmonic oscilla- 
tors. Because of poor screening the hole-apical inter- 
action is purely Coulombic, g^i'm. — n) = Ka/jm — n|^, 
a = x,y, z. To account for the experimental fact that 
z-polarized phonons couple to the holes stronger than 
the others [^ we choose ^ Ky = k^/V^. The direct 



hole-hole repulsion is Vc{n — n') 



so that the re- 



|n-n'| 

pulsion between two holes in the NN configuration is Vc- 
We also include the bare nearest neighbor (NN) hopping 
Tnn, the next nearest neighbor (NNN) hopping across 
copper Tnnn and the NNN hopping between octahedra 
T' 

According to Eq.(||), the polaron shift is given by the 
lattice sum (after summation over polarizations): 



En 



31.15 K^uj, 



m — n r 



m — n 



(13) 



where the factor 2 accounts for the two layers of apical 
sites. [For reference, Cartesian coordinates are n = {ux + 
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FIG. 1. Four octahedra of the copper-oxygen perovskite 
layer. Holes reside on the in-plane oxygens but interact with 
apical oxygens. 



1/2, Uy + 1/2, 0), m = (m^r, m^, h); rix-, ny,mx,my being 
integers.] The polaron-polaron attraction is 



^pa(n-n') 



AUJK 



2 + ^ n') • (m — n) 



(14) 



Performing lattice summations for the NN, NNN, and 
NNN' configurations one finds Fpa = 1.23 i^p, 0.80 E^p, 
and 0.82 Ep, respectively. Substituting these results in 
Eqs.(||) and we obtain the full inter-polaron inter- 
action: vj\iM — Vc ~ 1.23 Ep, vmnn ~ ~ 0.80 
^'nnn ~ ~ ^'^^ ^P' ^^'^ mass renormalization ex- 
ponents: G^^ = 0.38{Ep/Lu), G%j^pf = 0.60{Ep/uj) and 
G'nnn = 0.59{Ep/u;). 

Let us now discuss different regimes of the model. 
At > 1.23 £^p, no bipolarons are formed and 

the systems is a polaronic Fermi liquid. The po- 
larons tunnel in the square lattice with NN hopping 
t = Tnngxp{-0.38Ep/lu) and NNN hopping t' = 
TNNNe'>ip{-0.60Ep/uj). [Since G^^^ « G'^nn one can 
neglect the difference between NNN hoppings within and 
between the octahedra.] The single polaron spectrum is 
therefore 



i?i(k) ^ -Ep — 2t'[coskx + cosfcj^] 

± At cos{kx/2) cos{ky/2). 



(15) 



The polaron mass is m* = l/{t + 2t'). Since in general 
t > t' , the mass is mostly determined by the NN hopping 
amplitude t. 

If < 1.23 £^p, then intersite NN bipolarons form. 
The bipolarons tunnel in the plane via four resonating 
(degenerate) configurations A, B, G, and D, see Figure||. 
In the first order in -ffpcrt one should retain only these 
lowest energy configurations and discard all the processes 



FIG. 2. Top view on the perovskite layer. The apical sites 
are not shown. The four bipolaron configurations A, B, C, 
and D all have the same energy. Some possible single-polaron 
hoppings t' are indicated by arrows. Note that the bipolaron 
movement is first order in t' . 



that involve configurations with higher energies. The re- 
sult of such a projection is the bipolaronic Hamiltonian 

Hb = {Vc - 3.23 Ep) Y^iM^i + bIBi + G/Gi + d\Di] 
1 

- t'Y}AB\ + bICi + G/ A + dIAi + h.c] 
1 

n 

+ ClVxA + A-y^l+h.C.], (16) 

where 1 numbers octahedra rather than individual sites, 
X = (Ij 0)j and y = (0, 1). A Fourier transformation and 
diagonalization of a 4 x 4 matrix yields the bipolaron 
spectrum: 

£;2(k) = Vc- 3.23Ep ± 2t'[cos(fc^/2) ± cos{ky/2)]. (17) 

There are four bipolaronic subbands combined in a band 
of width 8t' . The effective mass of the lowest band 
is m** = 2/t'. The bipolaron binding energy is A = 
2£:i(0) - E2{0) = 1.23Ep -Vc-8t- At'. 

We have to emphasize that the bipolaron moves al- 
ready in the first order in polaron hopping. This remark- 
able property is entirely due to the strong on-site re- 
pulsion and long-range electron-phonon interaction that 
leads to a non-trivial connectivity of the lattice. This 
situation is unlike all other models studied previously. 
[Usually the bipolaron moves only in the second order 
in polaron hopping and therefore is very heavy.] In 
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our model, this fact combines with a weak renormal- 
ization of t' yielding a superlight bipolaron with mass 
TO** cx exp(0.60 E'p/cj). We recall that in the Holstein 
model TO** cx exp(2i?p/(jj). Thus the mass of the Frohlich 
bipolaron scales approximately as the cubic root of that 
of the Holstein one. 

At even stronger el-ph interaction, Vc < 1.16i?p, NNN 
bipolarons become stable. More importantly, holes can 
now form 3- and 4-particle clusters. Such clusters do not 
have resonant states and remain immobile in the first 
order in polaron hopping. The system quickly becomes 
a charge segregated insulator. 

The superconductivity window that we have found, 
l.\&Ep < Vc < 1.23i?p, is quite narrow. This indicates 
that the superconducting state in such systems is a subtle 
phenomenon which requires a fine balance between elec- 
tronic and ionic interactions. Too strong el-ph interac- 
tion leads to clustering, while too weak interaction cannot 
bind the carriers and the superconductivity is at best of 
BCS type. These considerations may provide additional 
insight into the uniqueness of one particular structure, 
the copper-oxygen perovskite layer, to HTSC. It also fol- 
lows from our model that superconductivity should be 
very sensitive to any external factor that affects the bal- 
ance between Vc and Ep. For instance, pressure changes 
the octahedra geometry and hence Ep and Vpa- Chemi- 
cal doping enhances internal screening and consequently 
reduces Ep. 

We now assume that the superconductivity condi- 
tion is satisfied and show that our 'Frohlich-Coulomb' 
model possesses many key properties of the underdoped 
cuprates. The bipolaron binding energy A should man- 
ifest itself as a normal state pseudogap with size of ap- 
proximately half of A Such a pseudogap is in- 
deed observed in many cuprates. There should be a 
strong isotope effect on the (bi)polaron mass because 
t^t' cx exp(— const -y/M). Therefore the replacement of 
O""^^ by O^^ increases the carrier mass Such an ef- 
fect has been observed in the London penetration depth 
of the isotope-substituted samples [0. The mass isotope 
exponent, am ~ dlnTO**/(ilnM, was found to be as large 
as a„i — 0.8 in Lai.895Sro.io5Cu04. Our theoretical ex- 
ponent is am = 0.3Ep/uj, so that the bipolaron mass 
enhancement factor is exp{0.6 Ep/uj) ~ 5 in this mate- 
rial. With the bare hopping integral T^nn = 0.2 eV we 
obtain the in-plane bipolaron mass to** ~ IOtoe. Cal- 
culated with this value the in-plane London penetration 
depth, Xab = [TO**/87rne^]^/^ ~ 316 nm (n the hole den- 
sity) agrees well with the measured one Xab — 320 nm. 
Taking into account the c-axis tunneling of bipolarons, 
the critical temperature of their Bose-Einstein condensa- 
tion can be expressed in terms of the experimentally mea- 
sured in-plane and c-axis penetration depths, and the in- 
plane Hall constant Rh as Tc = 1.64/ • {cRh / Xth>^c)^^^ ■ 
Here / « 1 and Tc, cRh, and A are measured in K, 
cm"^ and cm, respectively [E4|. Using the experimental 



Xab = 320 nm, Ac = 4160 nm, and = 4 x lO'^ 
cm'^/C (just above Tc) one obtains Tc = 31 K in striking 
agreement with the experimental value Tc = 30 K. The 
recent observation of the normal state diamagnetism in 
La2_2;Srj,Cu04 | p5[| also confirms the prediction of the 
bipolaron theory |26|. Many other features of the bipo- 



laronic (super)conductor, e.g., the unusual upper critical 
field, electronic specific heat, optical and tunneling spec- 
tra match those of the cuprates (for a recent review, see 
Ref. (13). ^ 

In conclusion, we have studied a model with strong 
long-range electron-phonon and Coulomb interactions. 
The model shows a reach phase diagram depending on 
the ratio of the inter-site Coulomb repulsion and the po- 
laronic (Franc-Condon) level shift. The ground state is 
a polaronic Fermi (or Luttinger) liquid at large Coulomb 
repulsions, a bipolaronic high-temperature superconduc- 
tor at intermediate Coulomb repulsions, and a charge- 
segregated insulator at weak repulsion. In the supercon- 
ducting phase, inter-site bipolarons are remarkably light 
leading to a high critical temperature. The model de- 
scribes many properties of the superconducting cuprates. 
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